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Abstract 

The circular orthogonal and circular symplectic ensembles are mapped onto free, non- 
hcrmitian fermion systems. As an illustration, the two-level form factors are calculated. 
February 7, 2008 

1 Introduction 

Dyson introduced the orthogonal, unitary, and symplectic ensembles of random matrices; differ- 
ent ensembles correspond to different physical problems according to the presence or absence of 
time reversal symmetry. The joint probability distribution function for the eigenvalues of matri- 
ces in these ensembles is equivalent to the statistical weight of a configuration of charges which 
repel each other with a logarithmic interaction at an appropriate temperature. By integrating 
over all levels except for some finite number, level-level correlation functions can be found[|^. 

Although the integrals can be carried out directly, it is of interest to find a physical 
system which reproduces the desired probability distribution function. This has previously 
been done using the Calogero-Sutherland modelfH, a model of interacting fermions. In this 
paper, physical systems will be presented which use non-interacting, but non-hermitian, fermion 
systems to produce the desired results. 
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2 Circular Ensembles 



The circular orthogonal and symplectic ensembles are defined by considering a system of N 
points on the unit circle in the complex plane Each configuration of points is weighted by a 
factor n 1 6*^' ~ e*^j |^, where /? = 1 for the orthogonal ensemble and /3 = 4 for the symplectic 

i<j 

ensemble. The points are located at the positions e*^*, for i = 1 to N. Both ensembles will be 
mapped onto fermion systems. 

The partition function will be given for arbitrary chemical potentials, which will permit 
the calculation of correlation functions, reproducing the Dyson formula Although this for- 
mula has been obtained previously by other means, and none of the final results are new, the 
fermion mapping for this problem provides a simple trick for calculating correlation functions. 



3 Orthogonal Ensemble 

First, we will consider the orthogonal ensemble. We wish to calculate the following partition 
function 

TV 

U J d0ill\''''' - e''^ (1) 

1=1 i<j 

Then, we will introduce a chemical potential and take functional derivatives to calculate corre- 
lation functions. 

Up to a sign, for given 9i, we may write 

We may write eq. (Q), up to a constant multiplier in front, as a correlation function in a 
two-dimensional free Euclidean fermionic field theory. We take the x coordinate to be the same as 
the 9 coordinate, and thus periodic with period 27r, and take correlation functions to be analytic 
in T + ix. We introduce the destruction operator ^(a;,r) which destroys a particle at (x,r), and 
the creation operator tp^XjT). The field ip has the action S = J dx drip^ {x,T){dr + idx)ip{x,T). 

Then we write eq. @ as 

{{f[ 0)){^t(o, -oo)}^/2{V.t(0, oo)}^/2) (3) 
1=1 
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where the operator {-0^(0,-0x3)}^/^ represents N/2 creation operators, shghtly separated from 
each other in space, and located very far away in time from the hne r = 0. The constant 
multipher relating eqs. (|2|) and (^) will depend on the precise separation between the creation 
points, and the precise distance between those points and r = 0. 

We may then forget about the two-dimensional field theory and write eq. as the 
expectation value 

N 

{y-\X{i^{x^)\V^) (4) 

1=1 

where the state |^^) is the state obtained by taking the vacuum state and adding particles to 
the N/2 lowest unoccupied energy states, and the state \V~) is the state obtained by taking 
the vacuum state and removing particles from the N/2 highest occupied energy states. This 
expectation value is the correct one since these two states are the ones projected out by the long 
time separation used in eq. (^). It should be noted that the constant factor relating eq. @ to 
eq. (H) is A/^-dependent, and, unless we take care to determine the constant, we will only be able 
to calculate correlation functions, and not the total free energy. 

The sign needed to relate eq. @ to the product in eq. (||) may be obtained by introducing 
an extra field r]. This field is x-dependent but confined to the line r = 0. For the ?7-field, x plays 
the role of time. This field has the action S = J dxT]rjx. Then the sign is simply given by 

N 

iUMxi) +rJix^)}) (5) 

i=l 

Then we write introduce a chemical potential f{x), weighting each configuration of points 

N 

by n /(^j)) s-iid write the partition function as 

i=l 

[dr]]{V~\e-f 
and integrate out rj to obtain 

{V~\e^ d-^d.x'f(x)tl;{x)f{x')ip(x')sign{x-x')^Y+\^ l^js^ 

where sign(x — x') is equal to 

J2e'^^''-'='^^ (8) 

k 

with k equal to half an odd integer. 
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It is convenient to rewrite eq. (^) for /(x) = 1 as 

{V-le" \V+) (9) 

This is clearly equal to 

(7V-l)/2 

(iV/2)! n Vk (10) 

k=l/2 

We may calculate the two-point correlation function of levels by taking two derivatives 
of eq. (|7[) with respect to f{x). The result splits into two pieces, eqs. ( pi] ) and (13). 

Calculating the correlation function of levels at points x and x' , one piece of the corre- 
lation function is obtained by inserting the operator 

-^(x)^(x')sign(x — x') (11) 

into the expectation value of eq. Expanding the ip operators in Fourier components 

a{k), a{k'), and noting that the result is non-vanishing only ii k = —k', the correlation function 
is 

- fce^'^C— ') ^ ie^K— ') (12) 
^ fe=_(Ar_l)/2 I 

where the sum is carried out over all half odd values of I. The factor of k arises due to a missing 
factor of ^ in the product of eq. (|lO|). 

The second piece of the correlation function is obtained by inserting the operator 

^ J dyd2/V(a^)^(y)sign(x - y) -0 (x') -0(7') sign (x - y') (13) 

We again expand each ip operator in a{k) operators, and then pair off operators a{k) with each 
other, so that two operators in a pair have opposite momenta. If we pair off the first with the 
second and the third with the fourth, the result is just the product of one-level correlators. 

If we pair off the first operator with the fourth operator, and the second with the third, 
we obtain 

{N-l)/2 {N-l)/2 

- Y (14) 

^ k=-{N-l)/2l=-(N~l)/2 
If we pair off the first operator with the third operator, and the second with the fourth, 
we obtain 

(7V-l)/2 (N-l)/2 
^ k=-(N-l)/2 l=-{N-l)/2 
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The average spacing between particles is 2tt/N . We rescale the x-coordinate to Nx and 
rescale the momenta k and / to k/N and l/N, to calculate the two-level correlator in units of 
average particle spacing. As becomes large, we may replace the sums by integrals, reproducing 
the Dyson formula. We calculate the m-th Fourier component of the form factor of the two-level 
correlator, where m = k + 1. 

Then, from equation (12) we obtain 



, , , .2\m\ + 1, 
1 + \m\ i-'o.y—^ — j 7) 



(16) 



'2\m\ - 1- 

From eqs. (14 Jl5|) we obtain a non-vanishing result only if \m\ < 1, in which case eq. (14) gives 



1 



and eq. (15) gives 



1 — \m\ + \m\ ln(2|m| — 1) 
So the form factor is given by eq. (p^) for \m\ > 1 and, for \m\ < 1, is equal to 



1 — 2\m\ + \m\ ln(2|m| + 1) 



(17) 



(18) 



(19) 



It is also possible to obtain the correlation functions of levels which belong to the same 
alternate series, by introducing separate chemical potentials for operators iprj and iprj. This may 
help relate this technique to the more standard method of integration over alternate variables 



4 Symplectic Ensemble 

For the symplectic ensemble, the desired partition function is 

N 

II deiUle'^^^-e'^^l^ (20) 

i=l i<j 

For given 0i, we may again write the weight of the configuration as a correlation function in a 
two-dimensional field theory. Here, the desired correlation function is 

N 

all V^(x„ 0)2){V;t(0, -oo)}^{V^t(o, 00)}^) (21) 

i=l 
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This is equal to 



which is 



where now N particles must be created at plus and minus infinite time instead of N/2 particles. 
The operator 7p(xi,0)'^ represents the operator 

lim-3-V(2;i + e,0)V(xi -e,0) (22) 
e-+o 2ie 

lim ^ e'^^+^'^''e'^^-''''^'a{k)a{k') (23) 
k,k' 

Y^e'^^+^'>^^^a{k)a{k') (24) 
k,k' ^ 

Taking \ V^) to be the state with particles added to the vacuum, and \V~) to be the 
state with N particles removed from the vacuum, eq. ( pOD becomes, up A^-dependent factors, 

(F^|e/'^^-^(^)'^(^)'|F+) (25) 

where a chemical potential f{x) has been introduced. For f{x) = 1, this is equal to 

2W-1/2 

A^! [| 2k (26) 

fc=l/2 

The calculation of the two-level form factor at this point is almost identical to that in the 
orthogonal ensemble. We must introduce the operator ^'0(x)^^/'(x')^ into the expectation value, 
expand in Fourier components, and pair off the annihilation operators in all possible ways. 

The non-trivial pairings of annihilation operators involve pairing the first with the third 
or the first with the fourth. These two possibilities give the same contribution. This yields 

^ i{k+l)(x-x') (fe - 0^ /27>, 
kl ^ 

which is zero for |m| > 2 and equal to, in the large N limit, 

\'ffi\ 

2 - \m\ + ^-^ln{\m\ - I) (28) 



for \m\ < 2. The factor of 1/4 in eq. (|2^ ) arises from the the factor of 1/2 in the operator 
inserted into the expectation value, times the factor of 2 from the two different pairings, times 
the factor of 1/4 from the two terms missing in eq. (|2^). 

A fermion mapping has been given for the orthogonal and symplectic ensembles. Al- 
though closely related to the technique of integrating over alternate variables, this technique 
provides a physical system, the fermionic field theory, which handles the various integrals and 
determinants. 
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